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ABSTRACT 

We examine the effects of the linear and the cosmological redshift-space distortions 
on the power spectrum of cosmological objects on a light cone. We develop theoretical 
formulae for the power spectrum in linear theory of density perturbations in a rigorous 
manner starting from first principle corresponding to Fourier analysis. Approximate 
formulae, which are useful properly to incorporate the redshift-space distortion effects 
into the power spectrum are derived, and the validity is examined. Applying our 
formulae to galaxy and quasar samples which roughly match the SDSS survey, we will 
show how the redshift-space distortions distort the power spectrum on the light cone 
quantitatively. 

Subject headings: cosmology: theory - dark matter - large-scale structure of universe 
- galaxies: distances and redshifts - quasars: general 

1. INTRODUCTION 

The wide-field redshift surveys of galaxies and quasars, like the Two-degree Field (2dF) and 
the Sloan Digital Sky Survey (SDSS), are key projects in the area of observational cosmology. 
These surveys provide three dimensional maps of the distribution of cosmological objects extended 
to high redshifts. The spatial distribution of the objects depends on various properties of the 
universe such as, initial density fluctuations, formation process of the cosmological objects, bias 
mechanism, luminosity function, cosmological parameters, and dark matter. The comparison of 
the observational data with theoretical predictions will put serious constraints on the cosmological 
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models. Especially the precise measurement of the power spectrum from the galaxy surveys is a 
promising project for the cosmological parameters and the initial density fluctuations, as well as 
the measurement of the cosmic microwave background anisotropies. 

Two-point correlation function and the power spectrum are the most used tools to quantify 
the clustering of the cosmological objects. From a theoretical point of view, these statistical 
quantities that are well defined on the constant-time hypersurface have been used so far. However, 
cosmological observations are carried out only on a light cone. This fact requires statistical 
quantities which properly incorporate the light-cone effect. For the 2dF and SDSS surveys the 
light-cone effect becomes important because of the depth of the samples. This importance of 
the light-cone effect was discussed by several authors (Matarrese et al. 1997; Matsubara, Suto, 
& Szapudi 1997; Nakamura, Matsubara, & Suto 1998; de Laix & Starkman 1998; Moscardini 
et al. 1998). Recently, a rigorous framework to incorporate the light-cone effect into the two-point 
correlation function and the power spectrum has been developed (Yamamoto & Suto 1999, Papar 
I; Nishioka & Yamamoto 1999, Paper 1; Suto et al. 1999 ; Yamamoto, Nishioka, & Suto 1999, 
Paper IE). 

The other important observational effect is the redshift-space distortion. Because cosmological 
surveys are carried out in redshift space, the peculiar velocity of sources distorts the distribution 
of the cosmological objects. This effect has been extensively investigated in linear theory (Davis 
&; Peebles 1983; Kaiser 1987; Hamilton 1998; Szalay, Matsubara, & Landy 1998) and also in 
a nonlinear regime (Cole, Fisher, & Weinberg 1994; Suto, et al 1999; Magira, Jing, & Suto 
1999). It has been also pointed out that the geometrical effect originating from the expansion of 
the universe causes another apparent distortion in the distribution of cosmological objects, i.e., 
cosmological redshift-space distortion (Alcock & Paczynski 1979; Ryden 1995; Matsubara & Suto 
1996; Ballinger, Peacock, & Heavens 1996; Suto, et al 1999; Magira, Jing, & Suto 1999; Matsubara 
1999). Though the authors in these previous works have focused on determining the cosmological 
parameters by using the cosmological redshift-space distortion, we note that the cosmological 
redshift-space distortion gives rise to a problem in determining the power spectrum because the 
cosmological parameters of the universe have not been well established. 

In this paper we focus on the power spectrum on a light cone, taking into account the 
redshift-space distortion effects. We derive formulae for the power spectrum in a rigorous manner 
and examine how these effects distort the power spectrum measured in wide- and deep-field 
surveys. We adopt the following strategy. We first define an estimator for the power spectrum 
which can be computed from a sample of cosmological sources on a light cone, according to the 
conventional Fourier analysis. Then we calculate the ensemble average of the estimator, which we 
define the power spectrum on a light cone. Comparing this with the conventional power spectrum 
defined on a constant-time hypersurface in real space, we show how the distortion appears in the 
shape of the power spectrum on the light cone quantitatively, modeling galaxy and quasar samples 
which roughly match the SDSS samples. 
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In the present paper we work in linear theory, which enables us to derive formulae in a 
rigorous manner. However, this limits the validity of our results only to the large scales because 
the nonlinear effect affects the clustering features on the small scales. So the extension to include 
the nonlinear effect is needed, and a plausible extension is discussed in Paper HI. 

This paper is organized as follows: In §2 we define the statistical estimators of the power 
spectrum and the two-point correlation function on a light cone. We will find that a familiar 
conventional relation between the power spectrum and the two-point correlation function, which 
is defined on a constant-time hypersurface, holds even for those thus defined on the light cone. In 
§3 we consider the linear redshift-space distortion effect on the power spectrum. Some parts in 
this section have been discussed in Paper HI. However, the approximate formula to incorporate 
the linear redshift-space distortion, which was introduced in Paper HI in an intuitive manner, is 
derived in the rigorous manner in the present section. And the validity of the approximate formula 
is examined in detail assuming galaxy and quasar samples. In §4 we consider the cosmological 
redshift-space distortion effect. Then the rigorous and the approximate formulae for the power 
spectrum on the light cone are derived. Then we examine the validity of approximation, and 
discuss the theoretical implications. §5 is devoted to summary and conclusions. Throughout this 
paper we use the unit in which the light velocity c equals I. 

For definiteness, we adopt notations similar to those in Paper HI, and use a superscript LC to 
indicate quantities on the light cone explicitly. The power spectrum without the superscript P(k) 
denotes the CDM power spectrum defined on a constant-time hypersurface at present. And we use 
subscripts (superscripts) R and S to indicate quantities in real and redshift spaces, respectively. 

2. DEFINITION OF THE POWER SPECTRUM 

In this section we first define statistical estimators for the two-point correlation function and 
the power spectrum of cosmological objects on a light-cone hypersurface defined by an observer 
in redshift space. We define the statistical estimators for the two-point correlation function and 
the power spectrum, independently, which can be computed when a set of a survey catalog is 
obtained. Then we will show that a familiar fundamental relation holds between thus defined 
two-point correlation function and the power spectrum. 

In a redshift survey, a position of a cosmological object is specified by the redshift z and 
the direction (unit) vector 7. Cosmological distance to the object is evaluated from the redshift 
z through some converting formula. And a map of the cosmological objects is obtained. The 
statistical quantities, the power spectrum and the two-point correlation function, are computed 
from the map. By introducing the variable s to specify the cosmological distance (the radial 
coordinate) in redshift space, we write the number density field of the cosmological objects 
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in redshift space n s (s,7) and the mean number density n s (s). [] Following the conventional 
treatment of the power spectrum (e.g., Feldman, Kaiser, & Peacock 1994), we introduce the 
number density field 

„/ x « S 0,7) -n^ ya (sn) , n x 
F{s,f) = — i , (1) 



d 3 sh S (sf 



where n^ yn (s,f) is the number density field of objects for a synthetic catalog without structure 
which has the same mean number density h (s) as that of the cosmological objects n °{s,~f). To 
be specific for the synthetic catalog, if we write n (5,7) and n-s yn (s,7) as 

n s (s, 7 ) = n s (s)[l + A s ( S , 7 )], (2) 

<n(*,7) = n s (s)[l + A s s yn ( S , 7 )], (3) 

respectively, we write 

[A^ n ( Sl ,7i)A^ n (s 2 ,7 2 )) = (A s ( Sl ,7i)A^ yn (s 2 ,7 2 ))= 0. (4) 

When the density field F(s,~f) is given, one may compute the following two-point statistics, 

J d i s 1 j d i s 2 F{s ull )F{s 2 , l2 )5^{s l - a 2 - R), (5) 

where s\ = si7j and s 2 = s 2 7 2 and R = \R\, and R = R/R. Note that this definition of the 
estimator for the two-point correlation function is slightly different from that in the previous paper 
(Paper I; Paper E). However, we will show that the difference is not practical and that the same 
formula is derived in the limit of the distant observer approximation (equation |2^] ) . 

Here ^ obs (R) is a mathematical expression for the two-point correlation function computed 
from a conventional data processing for a set of data n s (s,7) (for an observer). As is well-known 
as the problem of the cosmic variance, we can only predict the ensemble average of the statistical 
estimator, where the ensemble average means to average the estimator over many universes for 
different observers. Assuming that the selection function does not depend on the direction 7, the 
ensemble average of the estimator is written as 

J d 3 Sl j d 3 s 2 (F(s 1 , ll )F(s 2 , l2 )5^(s 1 -s 2 - R)), (6) 

which we define as the two-point correlation function on a light cone in redshift space. 

Next we consider the power spectrum. When the density field F(s, 7) is given, one can 
compute the Fourier coefficient of equation (||) 

jks 



/ d 3 s[n s (s,~f)-n s (s,~y)]e 
T{k) = / d 3 sF(s,~f)e lk - s = J - j 

r f ^3 „~s/„n2] 2 



d 4 sn*(sY 



(7) 



1 We assume infinitely thin binwidth for the 'observed' data set. 
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and the power spectrum may be computed as follows (Feldman, Kaiser, & Peacock 1994), 

P {k) = J l^^ (fc)l ' (8) 

where k denotes the wavenumber vector, k = \k\, and k = k/k. Similar to the case of the 
two-point correlation function, P obs {k) models a conventional estimation of the power spectrum 
computed from a set of data F(s,-f) for an observer. The ensemble average of the statistical 
estimator is written as 

P s LC (fc) = (p° bs (k)) = J ^i(|^(fc)| 2 ), ( 9 ) 

which we define as the power spectrum on a light cone in redshift space (see also Paper IE). 

When we adopt the above definitions, we show the familiar fundamental relation which holds 
between the two-point correlation function and the power spectrum, as follows. With the use of 
the relation 

_ 82 - R) = j MLjk-fr-*-*), (10) 



equation (||) is written as 



4vr J (2vr) 3 

= 2^2 / dkk 2 P^ c (k)j (kR). (11) 

The inverse Fourier transformation yields, 

P S LC (A;) = 4vr | dRR 2 £ c (R)j (kR). (12) 

Equations ([ll]) and (12) are the familiar formulae which hold for the conventional two-point 
correlation function and the power spectrum defined on a constant-time hypersurface. 



3. LINEAR REDSHIFT-SPACE DISTORTION 

The distortion due to the linear peculiar velocity field of cosmological objects is called as the 
linear redshift-space distortion. In this section we consider the linear redshift-space distortion 
under the assumption that the cosmological model is exactly determined. An error of assuming 
the cosmological model causes another distortion, i.e., the cosmological redshift-space distortion, 
which will be investigated in the next section. 



3.1. Two-point correlation function on a light cone 



In this subsection we first calculate the two-point correlation function, taking into account 
the linear redshift-space distortion. Formulation is developed within linear theory of density 
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perturbations based on the CDM scenario. Because the definition of the two-point correlation 
function is different from that in the previous paper (Papers H), derivation of the formula is 
slightly different from the previous one. We briefly summarize the calculation for definiteness. 

Throughout the present paper, we consider the spatially-flat Friedmann-Lemaitre universe, 
whose line element is expressed as 

ds 2 = a 2 (7?) [-dif + d X 2 + X 2 dn 2 2) ] , (13) 

where r\ is the conformal time, a is the scale factor, x is the radial coordinate, and dQ, 2 2 ^ is the line 
element of the unit two-sphere. We normalize the scale factor to be unity at present, 0(770) = 1. 
Then the Friedmann equation is 

where Hq = 100/ikm/s/Mpc is the Hubble parameter, Qo(= 1 — is the density parameter, and 
the dot denotes ^-differentiation. 

Since we locate a fiducial observer at the origin of coordinates (77 = 770, x = 0)j a cosmological 
object at x an d rj on the light cone hypersurface of the observer satisfies the simple relation 
77 = 7/0 — X- We introduce the radial coordinate r to express the three dimensional (real) space 
on the light-cone hypersurface, in which the position of a source is specified by (r, 7), where 7 
is an unit vector along the line of sight. Essentially, r is equivalent to x, however, we use r for 
mathematical conveniences and to represent that quantities are defined on the light cone explicitly. 
Then we denote the metric of the three dimensional space on the light cone as, 

dsl c = dr 2 + r 2 dnf 2) , (15) 

and we write the number density of sources on the light cone as n R (r, 7) , which is simply related 
to the comoving number density of objects at a conformal time rj and at a position, (x,7) , 

n R (r,7) = n(rj,xn)\v-*vo-r,x->r- (16) 

By using the mean observed comoving number density 72.(77) at time 77 and the density fluctuation 
of luminous objects A(t?, x>7)> we write 

n(v,X,l) =Hv)[ 1 + M r 1,X,l)]- (17) 
Then equation (]l6| ) is rewritten as 

n R (r, 7 ) =n R (r)[l + A R (r, 7 )], (18) 

where we defined 

n R (r) = n(7?)|^ _ r , A R (r, 7 ) = A(7/,x,7)|^o-r- )X -r- (19) 
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Note that the mean observed number density 77.(77) is different from the mean number density of 
the objects at 7/ by a factor of a selection function. 

Now we consider the two-point correlation function (^) in redshift space. The relation between 
the redshift space and the real space is specified as follows, 



n 



(s , 7) s A dsdVL~f = n R (r, rfr^drdthy, 



ra syn( s )7)s 2 c?sd^7 = nf ya (r,j)r 2 drdUj, 



(20) 
(21) 



and 



s = r + 8r(r, 7), (22) 
represents the apparent shift in the comoving radial coordinate due 



where 8r(r, 7) in equation (| 
to peculiar velocity. Sr(r, 7) is of order of the velocity perturbations, and the explicit expression 
is summarized in Appendix. Equation (|2l] ) defines the number density of the synthetic catalog in 
real space 7t R (r, 7). 

With the use of the above equations we obtain the following expression from equation (0), 

- 1 'do 



An 



dririd£l*y / dr2r 2 dCl~y 



d 3 rh R (r) 2 

x([7i R (n, 7l ) -n s R yn (n,7 1 )][n R (r 2 ,7 2 ) -«s R yn(^,72)] 
xS^ 3 ' (t*i + Sri — r 2 — <5r 2 — -R) 

where r + Sr = (r + Sr)~f, and we have assumed 



(23) 



J d 3 sh s (s) 2 = J d 3 rh R {rf 



(24) 



Note that the number density field of the synthetic catalog in redshift space (5,7) has no 



structure. Then the number density field in real space n R n (r, 7), which is defined by equation 
(|2l|), has a structure. Writing the number density field as 



<„(^7)=n R (r)[l + A R n (r,7)], 



syn\ 



the density contrast A R n (r, 7) is expressed as (Hamilton 1998), 



A R (r, 7 ) = A s b yn (r, 7 ) + 



syn 



where 



synv 



d nir) 

1 — 

dr r 



Sr(r,j), 



dlnr 2 ^' 



ir n (r) 



d In 'i 



(25) 



(26) 



(27) 



and we used Af yn (r, 7) = A~ yn (s,7), which holds within the linear order of perturbations. In the 
leading order of the perturbative expansion, equation (^) is rewritten as 



d 3 rn R (r 



R 



4tt 



drir(d£lj i I <ir 2 r2<ii7-y 2 fi R (r 1 )n K (r 2 ) 
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x([A> 1)7l )-( — + 



n 



<fri(n.,7i)] 



x[A R (r 2 , 72 ) - ( #- + ^V^(r 2 ,72)])^ (3) (n - r 2 - H) , 



(28) 



where we used equations 



<9r 2 r 2 
©, and (§). 

After straightforward calculations based on the CDM model with the linear biasing (equation 
8|), we obtain the following expression (see Appendix), 



drr 2 ^^) 



-i 



— / /_ dridr 2 rir 2 n R (ri)n R (r 2 ) 



■^2 I 'dkk 2 P(k)f[ 

i=l 



b(k;r] - t^D^q - n) 



j (kR) + P(k;r) - r 2 )I(k; R;r 1 ,r 2 ) + /3(/c;r? - n)I(k; R; r 2 , n) 



+P(k;r] - ri)/3(fc;77o - r 2 )J(k;R;r 1 ,r 2/ 



(29) 



where S denotes the region |ri — r 2 \ < R < r\ + r 2 , -P(A;) is the CDM power spectrum at present, 
Di(rj) is the linear growth rate normalized to be unity at present, I (k; R; r\ , r 2 ) and J(k; R;r\,r 2 ) 
are defined by equations (|A24j) and ( A25 ), respectively, and P(k;rj) is defined as 

1 din Di(rj) 



b{k;rj) d\n.a{rj) 



(30) 



3.2. Distant Observer Approximation 

As is shown in Paper H, equation (f29|) reduces to a simple form by applying the plane-parallel, 
or distant observer, approximation. In the case R <C 2r max , where 2r max is the size of the survey 
volume, we can use the approximation 

dr\dr 2 ~ / dr\ / dx (31) 
S J J-R 

where x = r 2 — r\, and I(k; R; n, r 2 ) and J(k; R; ri,r 2 ) are approximated as 

m,R;n,r,) = (32) 
J(kM) _ 3 ^-^ + M^, (33) 

respectively. Integration over x yields, 

f drr 2 n R (r) 2 C[R,z(r)] SOUIce 

£ C (R) ^ J - r—^— 2 . (34) 

drr n (r) 
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where we defined 

fft*(r)] B 



1 

2^ 



d/t/c 2 P(/c) J Di(??o - r) 2 6(A:; rj - r) 



r) + -(3(k;r] 
5 



jo 



(35) 



This formula indicates that the two-point correlation function on a light cone is obtained 
by averaging the correlation function at each time defined on a constant-time hypersurface 
£[R, z(r)] SOUTCe by weighting n R (r) 2 . Note that this formula has been rigorously derived only in 
the framework of linear theory of density perturbations. 



3.3. Power spectrum on a light cone 

The power spectrum is obtained from the two-point correlation function as described in 
section 2. By substituting equation (pii|) into (^) we compute the power spectrum on the light 
cone Pg jC (/c). Some aspects of the power spectrum on the light cone have been discussed in Paper 
M. Although a different method was adopted to compute the power spectrum in Paper HI, we find 
that the approach employed in this paper is rather useful. We derive the simple expression for the 
power spectrum, which was used without detailed verification in Paper HI, and discuss the validity 
of approximation in subsection 3.5. 

The approximate formula for the power spectrum is easily obtained by substituting equation 
|) into (0), 

P^ c (k) ~ a(k)P(k), (36) 



where 



J drr 2 h R (r) 2 b(k;r]Q — r) 2 Di(7]Q — r) 2 


l + lp{k;Vo-r) + -(3(k;r) -r) 2 
3 5 


J 


f drr 2 h R (r) 2 



a(k) 



This equation justifies the extended formula (20) in Paper HI. 



(37) 



In the case that the bias does not depend on scales of density fluctuation or k, i.e., 
b(k;rj) = b(rj), a{k) becomes a constant a. Thus in this case the light-cone effect and the linear 
redshift-space distortion are described by the constant a. In the small scales of large wavenumber 
of k, the nonlinearity distorts the power spectrum in redshift space, and produces an additional 
scale-dependence in the power spectrum on the light-cone. Importance of the nonlinearity in small 
scales is extensively discussed in Paper M (see also Magira, Jing, & Suto 1999). 
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3.4. Model of samples 

In subsection 3.5 we apply the formulae developed in the previous subsections to galaxy and 
quasar samples, and discuss the validity of approximation by comparing the approximate formula 
with the exact formula using numerical calculations. For that purpose, we need to model galaxy 
and quasar samples and bias. For a galaxy luminosity function, we adopt a B-band luminosity 
function of the APM galaxies (Loveday et al. 1992) fitted to Schechter function 

4>{L)dL = <p(±X exp (-£\ d (L\ (38) 



L* V L* \L 



with 0* = 1.40 x 10~ 2 /i 3 Mpc~ 3 , v = -0.97, and M% = -19.50 + 51og 10 /i. Then the comoving 
number density of galaxies at z which are brighter than the limiting magnitude B\\ m is given by 

h(z, < B hm ) = / <t){L)dL = <t>*Y[v + 1, x{B hra , *)], (39) 

JLIBu^.z) 



where 



x(B h 



L(Bi im , z) 



d L (z) 



1 2 



l/r^Mpc 



10 2.2-0.4B lim) (4Q) 



L* 

and T[v, x] is the incomplete Gamma function. We consider the galaxy sample in the range 
of redshift < z < 0.2, and adopt the B-band limiting magnitude 19 to match the SDSS 
spectroscopic sample. 

We also assume quasar samples which roughly match the SDSS quasar survey, where we use 
a selection function for B-band magnitude limited samples on the basis of the luminosity function 
by Wallington & Narayan (1993) and adopt the B-band limiting magnitude 20 (see also Nakamura 
h Suto 1997; Paper I). And we assume the depth of the survey volume, z max = 5. 

As for the bias, we adopt the simple scale-independent model (Paper IE), which is a 
phenomenological extension of the Fry's bias model (1996), 

6(t 7 ) = 1 + t -^— (&o-l), (41) 
[Di(r))]P 

where bo = b(rjo) is the bias parameter at present, and p is the constant which controls 
time-evolution. Note that this model reduces to the case of no bias when the constant bo is unity 
and that the case p = 1 is equivalent to the model by Fry (1996). 

We consider the SCDM(standard cold dark matter) and the LCDM(Lambda cold dark matter) 
models, which have (Oq, h, a$) = (1.0,0.0,0.5,0.56) and (0.3,0.7,0.7,1.0), respectively. We 



2 The terminology, 'exact' formula, means the exact expression rigorously derived from definition Q or (^) within 
the linear theory of density perturbation and bias. And the 'approximate' formula means its version in the distant 
observer limit. In a realistic situation, the nonlinearity of density perturbations and complex features of bias cause 
additional deformation of the power spectrum. (See also Paper IH, in which some aspects of the nonlinear effects are 
discussed.) 



- 11 - 



assume the Harrison-Zel'dovich initial power spectrum. The sets of the cosmological parameters 
are chosen so as to reproduce the observed cluster abundance (Kitayama & Suto 1997). And we 
use the CDM transfer function (Bardeen et al.1986; Sugiyama 1995) with Vl^h? = 0.015 in this 
section. Then the CDM power spectrum at present time is written as 

P(k) = BkT(k) 2 , (42) 

where 

T(k) = ln[l + 2.34g] 

{) 2.34< ? [l + 3.89g + (16.1g) 2 + (5.46<?) 3 + (6.71<7) 4 ] 1 /4 1 ' 

and 

' 2 ^\ 2 * (44) 

2.7 J n hexp[-n b - ^/2hn b /n ]hMpc~ 1, 

where B is a normalization constant. 



3.5. Validity of approximation and implications 

To compute the exact power spectrum on a light-cone, we first perform integration of equation 
(p9|), and obtain the two-point correlation function £g C (i?). Then we compute the power spectrum 
by using equation (|l^). As an example to show typical behavior of the two-point correlation 
function, Figure [l] plots the galaxy two-point correlation function in the LCDM model. Here we 
adopt the case of no bias, i.e., bo = 1 in (|4ll). The solid line, which represents the exact two-point 
correlation function (p9|), is compared with the dotted line, which represents the approximate 
formula (j34|). The exact two-point correlation function shows a good correspondence with the 
approximate formula at small R. However the deviation becomes significant at large R and the 
exact two-point correlation function drops down to zero at the point, where R becomes two 
times of the depth of the survey volume (diameter of the survey volume), i.e., 2r max (z max ). This 
is traced back to our definition of the two-point correlation function. Namely number density 
outside the survey region is zero, then a product of the number density of the two points being 
R > 2r max (z max ) must be zero. 

Now we discuss the power spectrum. Figure |2| plots the galaxy power spectra to show 
characteristic behaviors, in which we adopted the LCDM model and the case of no bias. The 
solid line, which represents the exact power spectrum, is compared with the dotted line, which 
represents the approximate formula (|36|). Comparing the exact and the approximate power 
spectra, we see that the approximation is fairly good for k 0.01/iMpc -1 . On the larger scales, 
the exact power spectrum approaches a constant value, and the deviation between the exact 
and the approximate formulae becomes large. This behavior originates from the fact that the 
corresponding wavelength becomes larger than depth of the survey volume at the large scales, 

r max , and that we cannot properly evaluate the power spectrum because of the finite size 
effect of the survey volume. For comparison, the CDM power spectrum defined on a constant-time 
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hypersurface at present is plotted (dashed line). Note that in this case the CDM power spectrum 
is identical to the galaxy power spectrum on the constant-time hypersurface because we have 
adopted the case of no bias. 

To quantify validity of the approximate formula, Figure || plots the power spectra for the 
galaxy and quasar samples, which are divided by the power spectrum on the constant-time 
hypersurface at present. The solid line expresses the exact power spectrum and the dotted line 
does the approximate one (|36| ) in redshift space. For comparison we show the case in real space, 
where we refer the 'real space' to the case f3(k; rf) = in the formulae. The short and the long 
dashed lines express the exact and the approximate power spectra, respectively, in real space. It 
is apparent from Figure |3| that the galaxy power spectrum is increased by the linear redshift-space 
distortion which dominantly contributes to the galaxy sample. The light-cone effect is not so 
effective for the galaxy sample because the sample is shallow. Nevertheless it can be notable that 
the light-cone effect decreases the amplitude of the power spectrum by order of several percent. 
On the other hand, the quasar power spectrum is significantly decreased because the light-cone 
effect is effective for such a deep observational field. However we should remind that the behavior 
of the power spectrum significantly depends on the time-dependence of the bias, which determines 
the amplitude of clustering at high-redshift (see also Figure |]). In each panel, we see that the 
approximation reproduces the exact power spectra fairly well for k ^ 0.01/iMpc^ 1 for the galaxies 
and quasars samples. Thus we conclude that equation (|^) provide a good approximate formula 
for the power spectrum on the light cone for k ^ 0.01/iMpc -1 . 

Here we briefly mention the nonlinear effects on the power spectrum. While the nonlinear 
evolution of density field enhances the amplitude, the fmger-of-God due to the random motion 
decreases the amplitude in redshift space. According to the result in Paper IE, the amplitude is 
decreased by order of several x 10% relative to its counterpart in linear theory for both galaxy 
and quasar samples, depending on the scale k and the cosmological model. This situation is 
clearly shown in Figure 2 in Paper IE, from which we see that the nonlinear effects are effective 
for k }t 0.1/iMpc -1 . Because we have not taken the nonlinear effects into account throughout the 
present paper, Figure || in the present paper does not show the decrease of the amplitude at the 
small scales. The difference between the Figure |j| and Figure 2 in Paper IH at small k is caused by 
the finite size effect of survey volume. 

Finally in this section we discuss dependence of cosmology. Figure |I] plots the factor a/b^ as 



a function of Q,q, where a is defined by equation fl37|). The solid and the dotted lines represent 
the factor ajb\ in redshift space and real space, respectively. Here we adopt the same model for 
the galaxy and quasar samples as in Figure ||. From Figure || it is apparent that the light-cone 
effect and the linear redshift-space distortion becomes more influential as S7o increases. This will 
be traced back to the linear growth rate of density perturbations. For the galaxy sample, the 
increase of amplitude due to the redshift-space distortion is more effective than the decrease due 
to the light-cone effect. In contrast to this, for the quasar sample, the light-cone effect and the 
time-evolution of the bias become important. 
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4. COSMOLOGICAL REDSHIFT-SPACE DISTORTION 

The crucial point of our investigation in the previous section is the assumption that we know 
the exactly correct cosmological model of the universe. Because the redshift z and the direction 7 
are the observable quantities of objects in the redshift survey, then we must assume a cosmological 
model to obtain a three-dimensional map on which the cosmological objects are plotted. Since the 
cosmological parameters of our universe have not been established completely, then we should be 
careful for this uncertainty. 

It is known that an apparent shape of distribution of cosmological objects is distorted by 
the cosmological redshift-space distortion. The cosmological redshift-space distortion has been 
discussed as a tool for the cosmological model (Alcock & Paczynski 1979; Ryden 1995; Ballinger, 
Peacock, & Heavens 1996; Matsubara & Suto 1996; Suto, et al. 1999; Magira, Jing, & Suto 1999). 
An incorrect assumption of the cosmological model affects proper estimation of the two-point 
correlation function and the power spectrum. In this section we focus on this point and we 
investigate how the error of the assumption of cosmological model affects the estimation of the 
power spectrum on a light cone. 



4.1. Basic Formulation 

In this section, we assume that the universe is the Friedmann-Lemaitre universe with a 
cosmological constant and that the 'correct' matter density parameter is Qq. We also assume that 
data processing is performed by assuming the Friedmann-Lemaitre universe with an 'incorrect' 
matter density parameter CIq. In this section we neglect the effect of the linear redshift-space 
distortion for simplicity. In this case, the real space is the universe with the 'correct' density 
parameter Qq, and is described by the same relations in the previous section. And we use the same 
notations from equation ( |i~3| ) to ( |l9| ) to describe the real space. On the other hand, the redshift 
space is the universe with the 'incorrect' density parameter CIq. To describe the redshift space, we 
write the line element 

ds 2 = d 2 (fj) [-dff + dx 2 + X 2 dn 2 2) ] , (45) 

where a{fj) is the scale factor normalized to be unity at present, i.e., a(fjo) = 1, fj is the conformal 
time, x is the radial coordinate, and and dfl 2 2 ^ denotes the line element of the unit two-sphere. 
Evolution of the scale factor a is specified by the similar equation as (p4|), 

where £l\ = 1 — Qq, and Hq = 100/ikm/s/Mpc is the Hubble parameter in redshift space. 
As in subsection 3.1, we assume that an observer is located at the origin of the coordinates 
(fj = Vo,X = 0)) an d the light-cone hypersurface of the observer satisfies the relation fj = fjo~ X- 
Then we introduce the radial coordinate s to denote the metric of the three dimensional 
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redshift-space on the light-cone hypersurface as 

«(2)- 



dsl c = ds 2 + s 2 dtt 2 (2) . (47) 



The relation between z and s is specified by 0(770 — s) = 1/(1 + z). When a set of data (2,7) 
for cosmological objects was obtained, we compute the number density field F(s, 7) defined by 
equation ((!]) with the number density field n s (s,7) and that for a synthetic catalog Ws yn (s, 7). 
And a map is constructed in redshift space (fl7[). Then we can compute the two-point correlation 
function and the power spectrum, according as the definitions © and @, respectively. We specify 
relations between the redshift space and the real space, as is done in subsection 3.1. The relation 
for the number density field is specified by the same equations as (^) and (21). The relation 
between s and r is specified by 

a( m - r) = 0(770 - s) = ——, (48) 

-L "1 Z 

where a{rf) and a(fj) are the solutions of equations ( |i~4| ) and fl46|), respectively. 

Now we calculate the two-point correlation function. With the use of equations, (|20|), (^3), 
and we find that (||) reduces to 

r d^t r r 

&W = i-i«J. 1 ; '- , (49) 

MSG) ^ 

where we used n {s)s 2 ds = n R (r)r 2 dr to write the denominator. 

After straightforward calculations similar to those in Appendix, we find that the above 
expression reduces to 



M j=1 



f s 2 _j_ s 2 _ £>2 \ 

x / dkk 2 P(k)b(k; r/ - ri)6(fc; r? - r 2 )j fcW rf + rf - nr 2 ^ 2 , ( 50 ) 



sis 2 

where M. denotes the region \s\ — S2I < < si + S21 and s\ and $2 are understood as functions 
of n and r 2 , respectively, which are specified by equation (|48"1) . Thus the formula ( |5C| ) is the 
exact two-point correlation function which takes the cosmological redshift-space distortion into 
account. Note that (|50|) reduces to the two-point correlation function in real space £^ C (-R), which 
is explicitly defined by ( |29| ) with setting j3{k\rf) = 0, when the 'correct' cosmological model is 
chosen as the redshift space, i.e., Q,q = £1$ and h = h. 
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4.2. Distant observer approximation 



We consider an approximate expression of ( |50[ ) in the simlar way to subsection 3.2. We 
introduce the variables y = s 2 — si and x = ri — r\, and set that y <C s\, S2 and x <C n.,?"2- In this 
case the argument of the spherical bessel function in (|5(]) is approximated as 



kJ r\ + r| — t'it^— ^— — — — - — k\ 



SlS 2 



1 



y 2 , 



where we used 



r 2 - r% dr\ 



y s 2 - si dsi ' 

Taking s\ and y as the integration variables instead of r\ and r 2 in (^), and we use the 
approximation , 



ds\ds2 



M 



R 



dsi / dy. 



(51) 
(52) 

(53) 



Introducing the new variable \i instead of y by fi = y/R, then (j5C|) is approximated as 



47r / dss z 



dr\' 
~ds) 



h K {rf 



i -l 



7T 



ds-\ 



dfj, 



dr\ x 2 x ^ 2 x 2 r 



si 



n R (r/ - ri) 1)1(770 - n) 



x y dkk 2 P{k)b{k\^-r x fiQ(km^, 



where we defined 



6 



+ 



(54) 
(55) 



Substituting equation (54) into (|l~2|), one obtains 



P s LC (fc) 



d^(^) (Ij) n R (r) 



x / dsis 2 



ds 



dsi 



™ R (% - n)Di(rio - n) 

2 







03 



(56) 



It might be more useful to rewrite this expression in terms of the mean number density per unit 
redshift and per unit solid angle dN(z)/dz, which is defined by dN(z) = h s (s)s 2 ds = n R (7')r 2 dr, 

r 1 2 



P S LC (£;) 



03 



where s is regarded as a function of 2 specified by the relation a(f]Q — s) = 1/(1 + z) 



(57) 
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4.3. Validity of approximation and implications 

We show how the cosmological redshift-space distortion distorts the power spectrum on a 
light cone. In this section we adopt the same galaxy and quasar samples as in subsection 3.4. We 
take f^o = 0.3 as the density parameter of the real space (the 'correct' universe). And we adopt 
CIq = 1.0 as that of the redshift space. In this section, we assume that h = h = 0.7 for the real and 
the redshift spaces for simplicity. f\ 

The exact power spectrum which takes the cosmological redshift-space distortion into account 
is computed from equation ( |i~2| ) by integrating the exact two-point correlation function (|50|). The 
approximate formula is given by equation ( |56| ) or (|57|). Figure ^ plots the galaxy power spectra. 
Figure || shows the case of the quasars. The solid and the dotted lines represent the exact and 
the approximate power spectra, respectively. It is apparent from the figures that the approximate 
formula shows the good correspondence with the exact formula for k ^ 0.01/iMpc -1 . The exact 
power spectrum approaches a constant value at the large scale, and the deviation of the exact 
formula from the approximate formula becomes large for k H 0.01 /iMpc" 1 . This is caused by the 
finite size effect of the survey volume, as mentioned in subsection 3.5. 

In Figures [| and ||, the short dashed line shows the power spectrum defined on the 
constant-time hypersurface at present P(k)b%. In lower panels, the power spectrum, divided 
by the power spectrum at present, Pg (k)/ P(/c)6q, is plotted. Although we here adopted the 
scale-independent bias model, the distortion of the power spectrum depends on the scale k. The 
location of the peak of the power spectrum is moved, due to the redshift-space distortion. Thus 
the cosmological redshift-space distortion distorts the power spectrum by shifting the wavenumber 
k, which is expected from the expression (|56|) or ([57]). This feature is in marked contrast to the 
case of the linear redshift-space distortion (see equation [^6| and Figure ||) . 

The long dashed line plots the case CIq = = 0.3, which is labeled by Pa (k). That is, the 
long dashed line shows the case when the cosmological redshift-space distortion is not effective, 
and it merely shows a contribution from the light-cone effect. Comparing the solid line and the 
long dashed line in the lower panels, we conclude that the cosmological redshift-space distortion is 
an important effect not only for the quasar sample but also for the galaxy sample in estimating the 
power spectrum. The cosmological redshift-space distortion affects the power spectrum at the 10 
percent level for the galaxy sample. For the quasar samples, the effect becomes more significant. 
Moreover we should note that the behavior significantly depends on the time-evolution of the bias 
model. 



3 The difference of the Hubble parameter can cause an additional difference in the power spectra which is not 
essential to the cosmological redshift-space distortion. Therefore we here choose h = h for definiteness. 
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5. CONCLUSIONS 

In the present paper, we investigated the power spectrum of cosmological objects on a 
light cone focusing on the redshift-space distortions in linear theory of density perturbations. 
We developed theoretical formulations to compute the power spectrum which properly takes 
into account the light-cone effect and the redshift-space distortion effects. The linear and the 
cosmological redshift-space distortions were considered in section 3 and 4, respectively. The 
formulae for the power spectra were derived in the rigorous manner starting from the first principle 
corresponding to the Fourier analysis. Using the distant observer approximation, we derived the 
approximate formulae (eqs.|}(| and f5(|), and examined the validity and the limitations. It was 
shown that the approximate expressions are useful to describe the effects properly. 

Applying our formulae to galaxy and quasar samples which roughly match the SDSS survey, 
we showed that the effects considered in the present paper become important for on-going wide and 
deep surveys like the SDSS and 2dF surveys. Let us summarize the results and the implications 
obtained in this paper. For the galaxy samples, the light-cone effect decreases the amplitude of the 
power spectrum several percentage points; the cosmological redshift-space distortion can distort 
the power spectrum by order of ten percent; and the linear redshift-space distortion increases 
the amplitude by order of several xlO percent, depending on the cosmological parameters. The 
power spectrum is rather insensitive to the time-evolution of the bias model for the shallow galaxy 
samples (Paper M). For the quasar samples, the linear redshift-space distortion can be a minor 
effect, and the light-cone effect and the cosmological redshift-space distortion are the influential 
effects. Of course the power spectrum significantly depends on the cosmological model and the 
bias model. 

In the present paper we have shown the validity of using the approximate formulas at the 
small scales. Our numerical calcuations have shown that the approximate formulas become 
invalid at the large scales for k ^ 0.01/iMpc -1 . Let us briefly discuss the physical meaning of the 
discrepancy between the exact formula and the approximate formula. The terminology, exact, 
means that the formulas are derived from the definitions of the statistical estimator for the 
two-point correlation function and the power spectrum, which corresponds to a conventional data 
processing. Then the exact formulas indicate how the finite size effect of the survey volume affects 
the two-point correlation function and the power spectrum. For the two-point correlation function 
the finite size effect decreases the amplitude of the correlation function at the large scales. For 
example the amplitude becomes zero for it! > 2r max , where 2r max is the diameter of the survey 
volume (see Figure 1). For the large separation (for large R), the sign of the two-point correlation 
function is negative, which implies the anti-correlation at the large scales. Thus the finite size 
effect decreases the anti-correlation at the large scales. We suppose that this causes the large 
amplification of the power spectra at the large scales (at small k) relative to its counterpart in 
linear theory. 

Because we worked in linear theory in the present paper, the behaviors on small scales (at 
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large k) are not correctly described, where the nonlinear effect causes an additional distortion 
of the power spectrum. The distortion of the power spectrum was discussed by extending 
the formula to incorporate the nonlinear effect in Paper HI, though the exact derivation was 
invalidated by the inclusion of the nonlinear effect in a strict sense. According to the result, the 
nonlinear effect causes the additional distortion of the power spectrum by order of several xlO 
percent at 0.1/iMpc^ 1 ^ k 1/iMpc -1 , depending on the cosmological model. In the present 
paper, the linear redshift-space distortion and the cosmological redshift-space distortion are 
discussed separately, however, investigations combining these effects are needed to present realistic 
theoretical predictions. (Suto, Magira, & Yamamoto 1999). Moreover we note that the bias 
model is the crucial problem to present precise theoretical predictions for the power spectrum 
(Moscardini, et al. 1998; Dekel & Lahav 1999; Taruya, Koyama, & Soda 1999; Tegmark & Peebles 
1998), though we here adopted a simple model. 
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APPENDIX 



A. Brief summary of the calculation £g C (i?) 

We here briefly outline the calculation of equation ( |29| ) . We use the almost same mathematical 
notations as those in Paper E, in which similar calculations have been given. The calcuations in 
Paper H will supplement arguments in this Appendix. We work within a framework of linear 
theory of density perturbations based on the CDM cosmological model. We expand the CDM 
density contrast 5 C and the velocity field v l c in terms of the scalar harmonics as follows (e.g., 
Kodama & Sasaki 1984) 

Uv,xn) = I" "ftEtotateT): (Ai) 

J ° l,m 

v t Mx,i)=J dkY, J ^r l yki m (x,i)K (A2) 

l,m 

where 34zm(x>7) is the normalized scalar harmonics 

3WX.7) =X l k ( X )Y lm (n 1 ), (A3) 

with 

Xi(x) = J^kjt(k X ), (A4) 

V 7T 

and Yi m ({l~y) and ji(x) are the spherical harmonics and the spherical Bessel function, respectively, 
and 3^m(x>7)'* denotes the covariant derivative of 34jm(x> 7) on the three-dimensional space. 
The Fourier coefficient S^lmi 1 !) satisfies 

4t+^l-|^41 = 0. (A5) 
In the Friedmann-Lemaitre universe, the growing mode solution is well known: 

4tfr)=4t(%)£>i(a), (A6) 

with 



p 'W = A V7 + 1 - n °j!,' fa ' U + ^(i-n.) J ' (A7) 

Here ^4 is a constant to be determined so that D\ is unity at present. 

We assume that the number density contrast of objects connects with the CDM density 
contrast by the scale-dependent bias factor b(k; rj) as, 

A klm ( V ) = b{k; V )8^ m ( V ). (A8) 
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Furthermore we assume that the velocity field of the objects, v(rj,x,y), is the same as that of 
CDM v c (rj,x,~f)- In this case the apparent shift of comoving distance of an object, Sr(r, 7), is 
related to the velocity as 



5r(r, 7) 



a(rj) 



H ^n + n A a(r/) 3 
With these assumptions we write 



7 • v e (r],x,y) 



where 



A R (r, 7 ) = [™ dkJ2d { kL(vo)b(k; V o ~ r)£>ifa, - r)y Wm (r, 7 ), 

roo . 

ffr(r, 7 ) = / dkJ26kL(Vo)f(Vo ~ r)D x { m - r)Ar 2 JW(r,7) l? 

J ° lm 



din a{rj) 

Substituting equations ( A10| ) and ( |A11|) into (|28|), we obtain 



(A9) 

(A10) 
(All) 

(A12) 



d 3 rn R (r) 2 



xn R (ri)ra R (r 2 )£>i(r?o - n)Di(r) - r 2 ) 



1 /• dft^ 

4"7T 



poo POO I \ 

/ *iE ^E(^U(%)42; m2 (%))^ imi (^7j^ m2 (%j 

J0 lim/ i 2 m 2 X ' 



n 

i=i 



b(ki;r] -n)-^ 2 V n 



where 



= f(rjo -r)— In 
ar 



r 2 h R (r)f(rj - r)£>i(r/ - r) 



- + /( o-r) — . 
<9r dr 2 



Using the relation 



i m=—l 



(A13) 

(A14) 
(A15) 



equation ( |A13| ) is written as 



d 3 rh R (r) 2 



1 r dn 



R 



47T 



drirfcZ07 i J dr2r\d£l~j ii 



xn K (ri)n K (r 2 )-Di(?7o - ri)£>i(% - r 2 ) 

c(c) 



limi " i 2 m2 



9 

x-n 

1=1 



&(A;i;?7o - n) - k 2 V r 



X7A3 / d3M7T E (-^JLdkrjYLMintfYZM^J 

x4vr ]T H) L2 jL 2 {kr2)Yl 2M2 {n k )Y L2M2 {Q l2 ) 
L 2 M 2 

x4vr ]T (i) L3 jL 3 (kR)YZ sM3 (n k )Y L3M3 (n A ) , 

L3M3 

where k = \k\ and k = k/k. Integration over ^7 2 , an d yields 

es°(«) = / d 3 rn R (rf ' J dnr? / dr 2 r|n R (n)n R (r 2 )Di(^ - n)£>ifa, - r 2 ) 



(A16) 



Zimi ' 



hm,2 



9 z 

x-n 

=1 

\2 



7T 



X 



(27T) 



^-3 y dkk 2 j h (kr 1 )ji 2 (kr 2 )jo(kR)5 hM 5^ 



mi,m2 • 



Because the Gaussian random fluctuations satisfy 



( r ?o)4 C 2 ) Z2m 2 ( r /o)) = - k 2 )Si 1 i 2 6 mim2 P(k 1 ), 



(A17) 



(A18) 



where -P(fc) is the CDM power spectrum at present, we obtain 



d 3 rh R (r) 2 



n Jo 



J dr x r\ J dr 2 r|n R (ri)n R (r 2 )Di(r/o - n)£>i (i]0 ~ r 2 ) 



— / dfcifc?P(A;i)X;(2i + l)Il 



b(h;vo -n) - k x 2 v r 



x y dkk 2 ji(k ri )ji(kr 2 )jo(kR) . 
By using the mathematical formulae, (Magnus et al. 1966) 



(A19) 



and 



{/ r 2_|_ r 2_^j2 \ 
ra^^^^J (|r 1 -r 2 |< J R<r 1 +r 2 ), ^ 
(i?< |n-r 2 |,i?>ri+r 2 ), 

^(2Z + l)P(cos 9)ji{kr 1 )j l (kr 2 ) = j (k^r 2 + r 2 - 2 n r 2 cos 0) , 



(A21) 



we find 



£ S LC (P)= [/ d 3 rh R (r) 

roc 2l. 

x y dkk 2 P(k)Y[ 



7rR J J s dridr 2 rir 2 n K {ri)n R {r 2 )Di(riQ -ri)Di(r/ - r 2 ) 



JO ^ V r l + r 2 ~ 2r l r 2 COS ( 



(A22) 
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where S denotes the region \n— r 2 \ < -R < ri-\-r 2 , and cos is replaced by cos = (r 2j i-r 2 — R 2 )/2r\r 2 
after operating the differentiations with respect to ri and r 2 . 

Using the formulae, (B17)-(B19) in Paper II, and omitting the second term in the derivative, 
i.e., approximating as T> r ~ f(j]o — r)d 2 /dr 2 , we finally have 



drr 2 n^"(r) 2 
1 



i -i 



1 



2vr 2 



2Rjj s 

2 



dridr 2 rir 2 n^'{ri)fi^'{r 2 ) 



dkk 2 P{k) Yl 



i=l 



jo(fc-R) + P(k;r) - r 2 )I(k; R;r 1 ,r 2 ) + /3(/c;r? - ri)I(k] R;r 2 ,n) 



+P(k;r]o - ri)P(k;rj - r 2 )J(k;R;n,r 2 ] 



(A23) 



where 



I(k;R;ri,r 2 ) 
J(k;R;ri,r 2 ) 



h(kR) j 2 (kR) {R 2 + r 2 
R 2 { 



J-,2 



j 2 (fei?) 

(kR) 2 



r\ + rf 



2r 2 
R 2 ^ 2 



1r\r 2 



+ 1 



+ 



j4(fcg) (R 2 + rl-rl\ 2 $R 2 + rl 



h(kR) 

kR 3 



2ri 
R 2 +r 2 



„2 \ 2 



2ri 

2 r>2 



+ 



2r 2 
« 2 + r 2 



R z + rf- ri R z + r 



r\ r\ + r\ 



2r 2 
-R 2 



r 2 



2r\r 2 



(A24) 



(A25) 



and (3{k;rjo — r) is defined by equation (30). Equation ( A23 ) is identical to equation (| 
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R [ h" 1 Mpc ] 

Fig. 1. — The two-point correlation functions on a light cone for galaxy sample in real and 
redshift spaces. The solid line represents the exact formula (29) and the dotted line represents the 
approximate formula (34). Upper pair curves correspond to the case that the linear redshift-space 
distortion is taken into account by using formulae (29) and (34). Lower pair curves correspond to 
the case of the real space by setting (3 = in (29) and (34). 
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Fig. 2. — The power spectra on a light cone for galaxy sample in real and redshift spaces. The solid 
line represents the exact power spectrum and the dotted line represents the approximate formula 
(36). Upper pair curves correspond to the case that the linear redshift-space distortion is taken 
into account, P^ c (k). Lower pair curves correspond to the case in real space, P^ c {k). (Here we 
refer the 'real space' to the case (3(k; rj) = 0. in the formulae.) The dashed line corresponds to the 
CDM power spectrum defined on a constant-time hypersurface at present P(k). 





6.0 




5.0 




4.0 




3.0 


eg 

o 




.Q 


2.0 






Q_ 


1.0 




0.0 


o 
_i 




Q_ 


1.2 




1.0 




0.8 




0.6 




0.4 




\SCDM galaxy (B<19z<0.2) - 
\ b =l.O 

A redshift-space exact formula 

\\ redshift-space approximate formula 

\\ ■ real-space exact formula 

- real-sapce approximate formula 




— I — I I I 1 1 1 1| 1 — I l i m i 1 — i i i 1 1 1 1 

SCDMQSO (B<20z<5.0) 
b =1.5 



■ ■ ■ ■ i i i i i i i i 



LCDM galaxy (B<19z<0.2) 



b =1.0 



"A 



-I — I I I III 1 — I I I I I 1 1| 1 — I I I I III 

LCDMQSO (B<20z<5.0) 
b =1.5 



■ ■ ■ i i i i i_ 



_i i 



10" 



10" 



10" 1 10"° 10" 

k [h Mpc" 1 ] 



10" 



Fig. 3. — Real and redshift space power spectra for galaxy and quasar samples, divided by power 
spectra of galaxy and quasar on a constant-time hypersurface at present, for SCDM and LCDM 
models. The solid line and dotted line represent the exact and approximate power spectra in redshift 
space, respectively. The short and long dashed lines express the exact and the approximate power 
spectra in real space, respectively. We adopted the bias model with bo = 1 for galaxies, and 
bo = 1.5, p = 1 for quasars. 
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Fig. 4. — The factor ajb\ for galaxy and quasar as a function of f^o- The solid and the dotted 
lines represent the factor in redshift and real spaces, respectively. The case of real space is defined 
by the case (5 = in the formulae. As for the bias, we adopt the model with 6o = 1 for galaxy, and 
bo = 1.5, and p = 1, 2 for quasar. 
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Fig. 5. — Galaxy power spectrum on a light cone which incorporates the cosmological redshift- 
space distortion. Here we adopted the galaxy sample in subsection 3.4, and (Qq = 0.3, h = 0.7), 
and (CIq = 1.0, h = 0.7). The case of no bias is considered. The solid and dotted line represent 
the exact and the approximate power spectra, respectively. The short dashed line plots the power 
spectrum on a constant-time hypersurface and the long dashed line labeled by P^ c (k) does the 
power spectrum assuming the correct cosmological model in redshift space, i.e., (f2o = = 0.3). 
Thus the long dashed line merely expresses the decrease of the power spectrum due to light-cone 
effect. The upper panel shows the power spectra, and the lower panel shows the power spectra 
divided by the power spectrum on the constant-time hypersurface at present, P(k). 
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Fig. 6. — Quasar power spectra on a light cone which take the cosmological redshift-space distortion 
into account. We adopted the quasar sample in subsection 3.4, and 0,q = 0.3, h = 0.7, and 
&o = 1.0, h = 0.7. Here the bias model with feo = 1.5, and p = 1, 2 is considered. The meanings 
of the lines are the same as Figure 5. 



